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We study nonclassical features in a number of spin-qubit systems including single, two and three qubit states, 
as well as an N qubit Dicke model and a spin-1 system, of importance in the fields of quantum optics and 
information. This is done by analyzing the behavior of the well known Wigner, P, and Q quasiprobability 
distributions on them. We also discuss the not so well known F function and specify its relation to the Wigner 
function. Here we provide a comprehensive analysis of quasiprobability distributions for spin-qubit systems 
under general open system effects, including both pure dephasing as well as dissipation. This makes it relevant 
from the perspective of experimental implementation. 
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I. INTRODUCTION 


A very useful concept in the analysis of the dynamics of classical systems is the notion of phase space. A straightforward 
extension of this to the realm of quantum mechanics is however foiled due to the uncertainty principle. Despite this, it is possible 
to construct quasiprobability distributions (QDs) for quantum mechanical systems in analogy with their classical counterparts 
M- These QDs are very useful in that they provide a quantum classical correspondence and facilitate the calculation of 
quantum mechanical averages in close analogy to classical phase space averages. Nevertheless, the QDs are not probability 
distributions as they can take negative values as well, a feature that could be used for the identification of quantumness in a 
system. 

The first such QD was developed by Wigner resulting in the epithet Wigner function (W) iHH. Another, very well known, 
QD is the P function whose development was a precursor to the evolution of the field of quantum optics. This was originally 
developed from the possibility of expressing any state of the radiation field in terms of a diagonal sum over coherent states 
IH3LH4I1 . The P function can become singular for quantum states, a feature that promoted the development of other QDs such as 
the Q function EHH as well as further highlighted the use of the W function which does not have this feature. These QDs are 
intimately related to the problem of operator orderings. Thus, the P and Q functions are related to the normal and antinormal 
orderings, respectively, while the W function is associated with symmetric operator ordering. It is quite clear that there can be 
other QDs, apart from the above three, depending upon the operator ordering. However, among all the possible QDs the above 
three QDs are the most widely studied. There exist several reasons behind the intense interest in these QDs. They can be used 
to identify the nonclassical (quantum) nature of a state & Specifically, nonpositive values of P function define a nonclassical 
state. Nonpositivity of P is a necessary and sufficient criterion for nonclassicality, but other QDs provide only sufficient criteria. 

A nonclassical state can be used to perform tasks that are classically impossible. This fact motivated many studies on non¬ 
classical states, for example, studies on squeezed, antibunched and entangled states. The interest on nonclassical states has been 
considerably amplified in the recent past after the advent of quantum information where several applications of nonclassical 
states, in particular, of entangled states, have been reported fl9T . Interestingly many of these applications have been designed 
using spin-qubit systems. 

Quantum optics deals with atom-field interactions. The atoms, in their simplest forms, are modeled as qubits (two-level 
systems). These are also of immense practical importance as they can be the effective realizations of Rydberg atoms 
Atomic systems are also studied in the context of the Dicke model ll22l[23h . a collection of two-level atoms; in atomic traps 
atomic interferometers |25ii], polarization optics |j2<3|, and have recently found applications in quantum computation ( I 27 I - I 32 I 


21 

24], 


and 


references therein) as well as in the generation of long-distance entanglement 13311 . All these would evoke the question whether 
one could have QDs for such atomic systems as well. Such questions, which are of relevance to the present work, would be 
closely tied to the problem of development of QDs for 5(7(2), spin-like (spin-j), systems. Such a development was made in 
& where a QD on the sphere, naturally related to the 5(7(2) dynamical group [35; 36], was obtained. There are by now a 
number of constructions of spin QDs B7l - l4lll . among others. 
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However, another approach, the one adapted here, is to make use of the connection of SU(2) geometry to that of a sphere. 
The spherical harmonics provide a natural basis for functions on the sphere. This, along with the general theory of multipole 
operators [42,43], can be made use of to construct QDs of spin (qubit) systems as functions of polar and azimuthal angles |44|. 
Other constructions, in the literature, of W functions for spin-1/2 systems can be found in j45[|46|] . among others. A concept 
that played an important role in the above developments, was the atomic coherent state [431, which lead to the definition of 
atomic P function in close analogy to their radiation field counterparts. Another related development, following ll48ll where joint 
probability distributions were obtained for spin-1 systems exposed to quadrupole fields, was a QD obtained from the Fourier 
inversion of the characteristic function of the corresponding probability mass function, using the Wigner-Weyl correspondence. 
This could be called the characteristic function or /' -function approach [49]. 

The fields of quantum optics and information have matured to the point where intense experimental investigations are being 
made. Both from the fundamental perspective as well as from the viewpoint of practical realizations, it is imperative to study the 
evolution of the system of interest taking into account the effect of its ambient environment. This is achieved systematically by 
using the formalism of Open Quantum Systems [f50l - l53h . 

In the present work, we investigate nonclassicality in a number of spin-qubit systems including single, two and three qubit 
states, as well as N qubit Dicke states and a spin-1 system, of importance in the fields of quantum optics and information. This 
is done by analyzing the behavior of the well known W, P, Q QDs on them. The significance of this is rooted to the phenomena 
of quantum state engineering, which involves the generation and manipulation of nonclassical states OH [55]]. In this context, it 
is imperative to have an understanding over quantum to classical transitions, under ambient conditions. Such an understanding is 
made possible by the present work, where investigations are done in the presence of open system effects, both purely dephasing 
(decoherence) l(56l 1531 , also known as quantum non-demolition (QND), as well as dissipation ll56l 15811 . These aspects of open 
system evolution have been realized in a series of beautiful experiments j59l |60f|. We also discuss the not so well known F 
function and specify its relation to the W function. Further, we expect this work to have an impact on tomography related issues, 
as borne out in [61], where a method for quantum state reconstruction of a system of spins or qubits was proposed using the 
Q function. Also, the Q function, studied here, can be turned to address fundamental issues such as complementarity between 
number and phase distributions Il62l - l64ll . under the influence of QND as well as dissipative interactions with their environment, 
as well as for phase dispersion in atomic systems [65!, .61]. Here, to the best of our knowledge, we provide, for the first time, a 
comprehensive analysis of QDs for spin-qubit systems under general open system effects. 

The plan of this paper is as follows. In the next section, we will briefly discuss the QDs that will be subsequently used in the 
rest of the work, i.e., the W, P, Q, and F functions. This will be followed by a study of open system QDs for single qubit states. 
Next, we take up the case of some interesting two and three qubit states as well as the well known N qubit Dicke model. We 
then discuss, briefly, QDs of a spin-1 system. These examples will provide an understanding of quantum to classical transitions 
as indicated by the various QDs, under general open system evolutions. Although QDs have been frequently used to identify 
the existence of nonclassical states [67], they do not directly provide any quantitative measure of the amount of nonclassicality. 
Keeping these in mind, several measures of nonclassicality have been proposed, but all of them are seen to suffer from some 
limitations lf68tl . A specific measure of nonclassicality is the nonclassical volume, which considers the doubled volume of the 
integrated negative part of the W function as a measure of nonclassicality l(69t] . In the penultimate section, we make a study of 
quantumness, in some of the systems considered in this work, by using nonclassical volume [69]. We then make our conclusions. 


II. DISTRIBUTION FUNCTIONS FOR SPIN (QUBIT) SYSTEMS 

Here, we briefly discuss the different QDs, i.e., the W, P, Q, and F functions, subsequently used in the paper. 


A. The Wigner function 

Exploiting the connection between spin-like, SU (2), systems and the sphere, a QD can be expressed as a function of the polar 
and azimuthal angles. This expanded over a complete basis set, a convenient one being the spherical harmonics, the W function 
for a single spin-/ state can be expressed as 0 

W (0, 0) = (2g!) 1/2 £ PkqYkq (0,0), (1) 

K.Q 

where K = 0,1,..., 2 j, and Q = —K, —K + 1,..., 0,..., K — 1, K, and 

= Tr{T' KQ p). 


Pkq 


( 2 ) 
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Here, Ykq are spherical harmonics and Tkq are multipole operators given by 


Tkq 


Z { 2 K+ i ) 1/2 

m,m' 


j K j 
—m Q m! 


I j, m){j , m' |, 


(3) 


where ^ ^ ^ ^ j = ( ^j+i —0'i TO iJ 2 W, 2 |j - to) is the Wigner 3j symbol (Tot] and (ji'rn\j 2 Til 2 1 j — to) is the Clebsh- 

Gordon coefficient. The multipole operators Tkq are orthogonal to each other and they form a complete set with property 
T^ = (—1)^ Tk,~q■ The W function is normalized as 


J W (0, </>) sin OdOcLcj) = 1, 

and W* ( 9 , <f>) = W (0 , </>). Similarly, the W function of a two particle system, each with spin-;/ is 

W (01, (/>!, 02, </>2) = (^rO Z £ PKiQiK 2 Q 2 Yk 1 Q 1 ($lj 0l) Yk 2 Q 2 (^2i ^ 2 ) : 

K!,QiK 2 ,Q 2 


(4) 


where 


Pk 1 q 1 k 2 q 2 = Tr {pT 1 KlQl Tj i2Q2 1 . Here, IF (6>i, 0i, 0 2 , <fe) is also normalized as 

W ( 0 i, </>i, 02, 02) sinsin0 2 (i0id(/)id02d(/ | 2 = 1 - 


/' 


Further, it is known that any arbitrary operator can be mapped into the W function or any other QD discussed here. In what 
follows, using the same notations we describe P, Q and F functions for single spin-/ state and for two spin-j particles. It may 
be noted that all the analytic expressions for the QDs given below are normalized. 


B. The P function 

In analogy with the P function for continuous variable systems, the P function for a single spin-/ state is defined as Hi 


and can be shown to be 


The P function for two spin-/ particles is 


P = J ded0P(e,4,) MM, 

(«,«= ew^ (»,«t a) 1/2 (-1 f- Q ( m -*as* +ii! ) 1/2 • 


P(01,01,02,02) = Z Z PK 1 Q 1 K 2 Q 2 Y Kl Q 1 (0i,4>i)Yk 2 Q 2 (02, fo) 

K u QiK 2 ,Q 2 

,, z' i\tri-Qi+tf 2 -Q2 l 1\ ( \/(2j-ifi)!(2j- if 2)!( 2 j+-f f i + 1 ) ! ( 2 J+-K'2 + l)! N \ 

X 1 Utt/ ^ (2j)!(2j)! ^ ■ 

Here |0, oi) is the atomic coherent state El and can be expressed in terms of the Wigner-Dicke states |j, m) as 

/ / c\ • \ 1/2 q r\ 

M = E ( m + jJ sin-' +m (-)cos^- m (-)e-^ + ^|j,TO). 


(5) 


( 6 ) 


(7) 


( 8 ) 


m=-j 


C. The Q function 


Similarly, the Q function for a single spin-/ state is 


Q 4 >) = </>)> 


(9) 
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and can be expressed as JH 

Q(M) = E PkqYkq ( 0 A) (^) 1/2 (~1)*- Q (2j +1) ( w J^ l K+ iy ) 


1/2 


K,Q 

Further, the normalized Q function for two particle system of spin-j 


(2j-K)\(2j+K+l)\ 

particles is 


Q = E E PK 1 Q 1 K 2 Q 2 Y Kl Q 1 (01,4>l) Yk 2 Q 2 (02, <fo) 

k 1 ,q 1 k 2 ,q 2 V > 


X (-1) 


Q1+K2 — Q2 


(2j)l(2j)l 


y /(2j-K 1 )\(2j-K 2 y.(2j+K 1 +l)\(2j+K 2 +iy. J ' 


( 10 ) 


( 11 ) 


D. The F function 


The F distribution function |49i] is defined using the relation between Fano statistical tensors and state multipole operators. 
Specifically, for a single spin-) state, it is defined as 


= Tj>kqYkq ( 0 . 0 ) ( A ) 1/2 £ ( ffi - gV )'" ■ 

Similarly, the normalized F function for a two particle, spin-j, Ei system is 

F $2) = E E Pk 1 q 1 k 2 q 2 Yk 1 q 1 {0i,(I)i)Y K2 q 2 {62,(I)2) 

k 1 ,q 1 k 2 ,q 2 


( 1 ^ 

( ( 2 j+K 1 +l)\( 2 j+K 2 +l)\ \ 1/2 

(477(2^1+^2) j 

v ( 2 j-K 1 y.( 2 j-K 2 y.{ju+i)} K i + K 2 J 


( 12 ) 


(13) 


To summarize, all the QDs discussed in this work are normalized to unity. They are also real functions as they correspond to 
probability density functions for classical states. The density matrix of a quantum state can be reconstructed from these QDs 
Q6D. One can also calculate the expectation value of an operator from them 14411. 

It would be appropriate here to make a brief comparison of the QDs, discussed above, with their continuous variable counter¬ 
parts. The coherent state and thereby the displacement operator D, which generate coherent states from vacuum and is usually 
expressed as D((f> ) = e'® 0,1 '® a with a, a' being the annihilation and creation operators of the given Fock space, respectively, 
play a central role in these considerations. Thus, for example, the Wigner function, associated with a state p is the symplectic 
Fourier transform of the mean value of D in the state p , leading to the standard representation of the Wigner function as the 
Fourier transform of the skewed matrix representation of p, 

1 f°° £ £ 

W(rj 1 rj)=—J d£(a~ ||p|a + |}e l?b , (14) 


with a and b being real and ?/ = (a + ib). Similarly, the P function is associated with the diagonal representation of the state 

p in terms of the coherent state, while the Q function is related to the expectation value of p, with respect to coherent states. 

The multipole operators Tkq © play a pivotal role in the construction of QDs of spin systems, discussed here. These 
operators are extensively used in the study of atomic and nuclear radiation and can be shown to have properties analogous 
to those of the coherent state dispalcement operator D for usual continuous variable bosonic systems |44j|. In this sense, the 
properties of spin QDs are analogous to those of their continuous variable counterparts, with the atomic coherent state playing 
the role of the usual coherent state. Thus, for example, even though, for the spin (qubit) systems, both P and W QDs are 
witnesses of quantum correlations, in the sense that their negative values indicate quantumnes in the system, it is possible for a 
scenario wherein the P function is negative and W is positive, but not viceversa. The Q function is always positive, while the F 
function is same as the W function for spin-j systems, as shown below. 

Before proceeding further, it is worth noting here that for all the spin-j states (qubits), single or multi-qubit, the W and F 
QDs are identical. Specifically, for the single qubit case, the W function is 

Wx(9,<f>) = -±=Yi M K’Q 

V 27r K,Q 


where Mk,q (0, <j>) = PkqYkq (0, <j>), while, the F function is 


Fi(0,<t>) = -±=Y l M K’Q( e ’ < fi( 


(2 + K)\ \ 1/2 


3 K (1 — K)\ 
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where the term inside the brackets with square root is 2 for both the values of K (i.e., 0 or 1). Similarly, for two spin-i states, 
the W and F functions are 

Fi 1 ( 6 > 1 , </)i, 02, 02) = ^ m k 1 q 1 k 2 q 2 (01,01, 6 * 2 ,02) 

k 1 ,q 1 k 2 ,q 2 

= 57 J2 M Kl Q 1 K 2 Q 2 (01,^1,02)02) 

k 1 ,q 1 k 2 ,q 2 

= W h i (0i,0i,0 2 ) 02), 

where M KiQiK2 q 2 (9 1 , 0i, 6 2 , fa) = PKiQiK 2 q 2 Yk 1 q 1 (0i, </>i) Y K2 q 2 ( 9 2 , fo ), and the term in the brackets is 4 for all the 

possible values of K i and K-i- This can further be extended for higher number of spin--j states. 

Since the W and F functions are the same for spin-i systems, we will not discuss the evolution of the F function of these 
systems. 


( (2+K 1 )\(2+K 2 )\ 

- 


1 + « 2 ( 1 -Jf 1 )!(l -K 2 ) 


0 


i/z 


III. DISTRIBUTION FUNCTIONS FOR SINGUE SPIN-1 STATES 


Here, we consider single spin-1 states, initially in an atomic coherent state, in the presence of two different noises, i.e., QND 
lf56L l57h . which are purely dephasing, and the dissipative SGAD (Squeezed Generalized Amplitude Damping) |56j, 58j| noises. 
For calculating the QDs, we will require multipole operators for j = 1 and m, m! = ±1, giving K = 0 and 1. For I\ = 0, 


<3 = 0, and for K = 1, Q = 1, 0, —1. Using these, the multipole operators Tkq can be obtained as Too 


l 

V2 


1 0 
0 1 ’ 


Tn 


0 0 

-1 0 ’ 


Tio 


i 

V2 


1 0 
0 -1 


and Ti_i 


0 1 
0 0 


A. Atomic coherent state in QND noise 


The master equation of the system interacting with a squeezed thermal bath and undergoing a QND evolution IH i 


is 


Pnm (t) = -JJ (E n - E m ) + if] ( t ) (El - Ell) - ( E n - E m ) 2 7 (t) p nm ( t ) , 


(15) 


where E n s are the eigenvalues of the system Hamiltonian in the system eigenbasis |n), which here would correspond to the 
Wigner-Dicke states |4I1 ; 


k 


and 

7 W = coll 1 ( '2 k ) I (e lWkt — l) cosh( r &) + _ D s j n h e 2l ^k\- _ 

k k 1 ' 

Here, f3 = jjpp* and k/j is the Boltzmann constant, while 77 and 4 > k are the squeezing parameters. The initial density matrix for 
the atomic coherent state is 


p(0) = \a, /3)(a, /3\, 

where \a, /3) is given by Eq. ©. The density matrix (1 1 6[ i in the presence of QND noise at time t becomes 

pjm.fr (t) = ( 0 ), 

where 


For j = 1, the initial density matrix is 


Pjm,jn{ 0) = {j, m\p (0) |j, n) 

= {j,m\a,P)(a,P\j,n). 


P( 0) 


sin 2 (f) \e 1/3 sin a 
le l/3 sina cos 2 (f) ’ 


(16) 

(17) 


(18) 


(19) 
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Figure 1: (Color online) Plot (a) shows the variation of all QDs with time (t) for single spin-| atomic coherent state in the presence of 
QND noise with bath parameters 70 = 0.1, oj c = 100, squeezing parameters r = 0, a = 0, and oj = 1.0 at temperature T = 1, and 
a = ^,/3 = 9 = (j) = j, in the units of h = ks = 1. The smooth (blue) line, dashed (red) line and dot-dashed (magenta) 

line correspond to the W, P and Q functions, respectively; (b) and (c) show the variation of W and P functions with time for different 
temperatures T = 0, 1 and 2 by smooth (blue) lines, dashed (red) lines and dot-dashed (magenta) lines, respectively. 


which in the presence of QND noise becomes 


P(t) 


Bin 2 (f) 




sin ct 


\ e -i“t e -(M 2 7(t) e -i/3 s in a 
cos 2 (f) 


( 20 ) 


Here, we consider the case of an Ohmic bath for which analytic expressions for 7 (£), both for zero and high temperatures, can 
be obtained fUzti- These are functions of the bath parameters 70 and lo c as well as squeezing parameters r and (!>, with (f> = auj 
and a is a constant dependent on the squeezed bath. Now, using multipole operators, mentioned above, analytic expressions of 
the different QDs can be obtained. For example, we have obtained the W function for a qubit, starting from an atomic coherent 
state, in the presence of QND noise as 

W(9,cj)) = ^1 — \/3 cos a cos 9 + ■y/Se cos (/? + ojt + <j>) sin ct sin 6^j , (21) 

while, the corresponding P and Q QDs are obtained as 

P{9 1 (j)) = ^1 + 3 cos ol cos 9 + 3e — cos (/3 + ojt + </>) sin ct sin 9^j , (22) 

and 

Q(9,cj>) = ^1 + cos a cos 0 + cos (/3 + cot + (f>) sin a sinl^ , (23) 

respectively. All the QDs calculated in Eqs. (|2T1 >- (|23| > can be used to get the corresponding noiseless QDs for the same system 
and this also serves as a nice consistency check of the calculations. The variation of the QDs, Eqs. (|2T1 i-(| 23|>. for some specific 
parameter values are shown in Fig. Q]a-c, where the effect of the presence of noise on the QDs can be easily observed. Both 
the P and W functions are found to exhibit negative values indicative of quantumness in the system. Also, in Fig. |T|b-c, we 
can see that with an increase in temperature T, the QDs tend to become less negative, which is an indicator of a move towards 
classicality, as expected. Interestingly, in Fig. [I] a, we do not observe any zero of the Q function which implies that the Q 
function does not show any signature of nonclassicality in this particular case. The oscillatory nature of the QDs for atomic 
coherent state when subjected to QND noise can be attributed to the purely dephasing effect of the QND interaction. That is, this 
process involves decoherence without any dissipation. At temperature T = 0 decoherence is minimal and hence an oscillatory 
pattern is observed in the depicted time scale. With increase in T, resulting in increase in the influence of decoherence, these 
oscillations gradually decrease. 


B. Atomic coherent state in SGAD noise 

Now, we take up a spin j = 2, starting from an atomic coherent state, given by Eq. ( | 1 9k evolving under a Squeezed Gener¬ 
alized Amplitude Damping (SGAD) channel, incorporating the effects of dissipation and bath squeezing and which includes the 
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well known amplitude damping (AD) and generalized amplitude damping (GAD) channels as special cases. The Kraus operators 
of the SGAD channel are [58] 


E 0 = 


o 

o 1 


Ei = y/p 


0 0 


Vm o 


e 2 = v 1 — p 


VH4) 

0 


V 1 -"(t). 


E 3 


V 1 -p 


0 


0 


(24) 


where A = ± {1 — (1 — p) [/x + v) - exp (- 7o (2JV + 1) t)} , p = 2 ^i-?) sinh^of^+i'U) exp ( 2iV + X ) *) > and 
zz = ( 1 _p)^^+i) {1 — exp (—70 {2N + 1) t)} . Here, for convenience we have omitted the time dependence in the argu¬ 
ment of different time dependent parameters (e.g., A (t), p(t), v{t), etc.) in the Kraus operators of SGAD noise. Here, 
70 is the spontaneous emission rate, a = sinh(2r) (2 Nth + 1 ), and N = N t h {cosh 2 (r) + sinh 2 (r)} + sinh 2 (r), with 
Nth = 1 / {exp ( huj/kBT ) — 1} being the Planck distribution. Here, r and the bath squeezing angle (£ (t)) are the bath squeez¬ 
ing parameters. The expression for p in the above equations has an analytic, though complicated, expression, and we refer the 
reader to |[58ll for details. Application of the above Kraus operators to the initial state results in 


Pit) = ,S^(t)p(0)Bj(t). 

i—0 


The state at time t can be obtained as 


Pit) 


P 11 P 21 

P 21 1 — P 11 


(25) 


where 


P 11 = \ {1 — P + v ~ P (A — p + v) + (—1 + n + v + p (A — p - v)) cos a} , 

and 

P 21 = | sin a {(1 - p) 'Jiwe-W+S + P Vl^Xe* + (1 - p) V(1 - Ai) (1 - i/)e i/S } . 

Using this density matrix, we can calculate the evolution of the different QDs, in a manner similar to the previous example of 
evolution under QND channel, leading to 

W [l + v / 3{— p + v — p (A — /i + v) + (— 1 + /i + v +p (A — p — v)) cos a} cos0 

+ v/3 (jpV 1 - A + (1 -p) y^l - p) (1 - zz)| cos (/3 + <j>) (26) 

+ (1 - p) y/jlv cos (P+ £-(/>)) sin a sin 0] ; 

P{6,(j>) = [1 — 3 {—p + v — p (A — p + v) + (—1 + p + v +p (A — p — v)) cos a} cos# 

+ 3 (jpVl — A + (1 — p) -y/(l — /x) (1 — t/)| cos (/3 + <j>) (27) 

+ (l-p)V^ c °s(/3 + ?-0)) sin a sin d\ ; 

and 

Q(Q,<I>) = -E [l - {-p + v - p (A - p + v) + (-1 + /i + v +p(A - p - u)) cos a} cos6» 

+ 1 — A + (1 — p) \J (1 - p)(l-v)\ cos (/? + <\>) 

+ (1 — p) y/Jwcos (/3 + £ — <j))) sinasind] . 


( 28 ) 
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Figure 2: (Color online) The variation of all the distribution functions with time (f) is shown together for a single spin -77 atomic coherent state 
in the presence of the SGAD noise for zero bath squeezing angle in units of h = Icb = 1, with u = 1.0, 70 = 0.05, and a = j, /3 = ^, 6 = 
(f> = In (a) the variation with time is shown for temperature T = 3.0 in the absence of squeezing parameter, i.e., r = 0. In (b) the effect 
of the change in squeezing parameter for same temperature, i.e., T = 3.0 is shown by using the squeezing parameter r = 1.0, keeping all the 
other values as same as that used in (a). Further, in (c) keeping r = 1.0 as in (b), the temperature is increased to T = 10 to show the effect of 
variation in T. In (c) time is varied only up to t = 5 to emphasize the effect of temperature. In all the three plots, smooth (blue), dashed (red) 
and dotted-dashed (magenta) lines correspond to the W, P and Q functions, respectively. 


The variation of all the QDs with time (t) for some specific values of the parameters is depicted in Fig. [2] which incorporates 
both temperature and squeezing. A comparison of the Figs. [2] a and[2]b brings out the effect of squeezing on the evolution of 
QDs. Further, it is easily observed that with the increase in T, the quantumness reduces. An important point to notice here, is that 
if we make the noise parameters zero, i.e., in the absence of noise, the different QDs given by Eqs. (I26l >- (f28] >. reduce to a form 
exactly equal to the corresponding noiseless QDs obtained for QND evolutions (Eqs. (|2TI >- (|23] >). Also, results for generalized 
amplitude damping channel can be obtained in the limit of vanishing squeezing, i.e., for p(t) = 0 and A (t) = v (t), while 
corresponding results for QDs under evolution of an amplitude damping channel can be obtained by further setting T = 0, and 
p = 1. Further, it would be apt to mention here that the oscillatory nature of the QDs for an atomic coherent state evolving under 
QND noise is not seen here. This is consistent with the fact that the SGAD noise is dissipative in nature, involving decoherence 
along with dissipation. 


IV. QDS FOR MULTIQUBIT SYSTEMS UNDERGOING QND AND DISSIPATIVE EVOLUTIONS 

Now, we wish to study the evolution of QDs for some interesting two and three qubit systems under general open system 
evolutions. We will also take up the well known TV-qubit Dicke model. In each case, we study the nonclassicality exhibited by 
the system under consideration. 


A. Two qubits in the presence of QND noise 


The density matrix for a system of two qubits in QND interaction with a squeezed thermal bath, as obtained in Ref. HHu 


IS 


P \in,3n}^> = eX P [*{©{!„,(0 W}] eX P 




(29) 


where p^ ? , (t) is the two-qubit reduced density matrix obtained by tracing out the bath (reservoir) degrees of freedom and 
has the matrix representation (£2, *i| p s it) |j2, Ji), and {i n ,jn} stands for In this model, the system-bath coupling 

is dependent upon the position of the qubit, resulting in the classification of the dynamics into two regimes: (a) Localized 
model, where the inter-qubit spacing is greater than or of the order of the length scale set by the bath, and (b) Collective model, 
where the qubits are close enough to experience the same bath lt7~ill . Here, for the sake of brevity, we will provide details of 
the localized model only. The terms 0{ 7 ;„, Jn } (t), A {i„ ,j n } (t) and r?? 3 j (t) have different expressions in the localized and 
collective models. The superscript sq indicates that the bath starts in a squeezed thermal initial state. For convenience, the two 
particle index in Eq. ( l29l ) is denoted by a single 4-level index in the following manner: 

——, —- = 0; = -,-- = 2- -,- = 3. 

2 2 2 2 2 2 2 2 
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All the sixteen terms, of the density matrix, can be analytically calculated for a given initial state p s (0). In the localized model, 
considered here, the density matrix is obtained using the symmetry of the density matrix p s (t), i.e., symmetries between the 
matrix elements and hermiticity of the density matrix, and the expressions of different terms in Eq. m El. The elements 

P32 (t) = P23 (t) ~ P01 (t) = P10 (t) i 


are obtained using 


©32 (t) = ©01 (f) = —©23 (t) = —©10 (t) 

= f Q dull (ui) S (ui,t) cos uit s , 

and 

A32 (t) = A01 (t) = — A23 (t) = — A10 (t) 

= — / 0 °° dull (w) C (w, t) sinwf s , 

where I (ui) is the bath spectral density. In the Ohmic case considered here, I (w) = where 70 and 0 o c are bath 

parameters. We have S (w,t) = ujt ~ s ) 2 nujt . C (u,t) = wi s = k.r mn , where r mn is the inter-qubit spacing and t s is 

the transit time introduced for the purpose of expressing the coupling of the system to its bath in the frequency domain. The 
diagonal elements of the density matrix are 


P s aa 0) = Paa (0), where a = 0,1, 2,3, 

implying an unchanging population, a characteristic of QND evolution. Also, 

P21 (t) = P12 (t) = P12 (t) f 
P 30 ft) = POS (^) = P 03 (t) > 

i.e., these elements are purely real and for these 0 (t ) = 0 = A (t). The remaining elements are 

P31 (^) = P13 (t) = P02 (t) = P20 (t) j 


and for their calculation we need 


©31 (t) = ©02 (f) = —©13 (f) — —©20 (t) 

= f Q dull (ui) S (ui, t) cosuitg, 


and 


A 3 I (t) = A 02 (f) = — A 13 (f) = — A 2 0 (t) 

= / Q °° dull (ui) C(u>,t) sin u>t s . 

For the determination of all these elements of the density matrix at time f, we also need T S9 (t) which have complex expressions 
and can be seen in 0 . Once the density matrix at time t is obtained, the corresponding QDs can be obtained from the 
prescription discussed above. However, getting analytic expressions for different QDs, here, is a cumbersome task; hence we 
will resort to numerically plotting them. Without loss of generality, we consider here our initial state to be such that all the 
sixteen elements of the density matrix at time t = 0 are 0.25. The effect of QND interaction on different QDs obtained for this 
particular choice of initial state is shown in Fig. [3] As can be seen from the figure, both the P and IT functions exhibit negative 
values, indicative of quantumness in the system, for some time after initiation of the evolution, before becoming positive, due to 
dephasing caused by the bath. As expected, the P function is a stronger indicator of quantumness than the W function, while 
the Q function is always positive, by construction. 


B. Two qubits under dissipative evolution 

Here, we study the evolution of QDs, for two qubit systems, undergoing dissipative evolution, first interacting with a vacuum 
bath, T = 0 and zero bath squeezing, and then under the influence of a squeezed thermal bath, finite T and bath squeezing. 
Here, we will make use of the results worked out in Ref. l(72h . 
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QDs 



Figure 3: (Color online) The plot shows the variation of all the QDs with time for two qubits undergoing QND evolution with 70 = 0.01, oj c = 
100, k.r a b = 0.05, squeezing parameters r = 0.05, a = 0, which implies zero squeezing angle, and temperature T = 2 with #1 = -|, d 2 = 
= 7 t, 02 = 7 in units of h = fes = 1. The smooth (blue), large dashed (red), and small dashed (magenta) lines correspond to W, P 
and Q functions, respectively. 


1. Vacuum bath 

The density matrix, in the dressed state basis, can be used for calculating different QDs. We consider the initial state with one 
qubit in the excited state \e\) and the other in the ground state \g 2 ), i.e., \e\) | g 2 ). The two-qubit reduced density matrix is given 
by 


Pit) 


Pee {t) pes (t) 
P*es(t) Pss{t) 

Plait) P * sa (t) 

l Pig it) Pig it) 


Pea it) Peg it) 

Psa it) Psg it) 

Paa it) pag (t) ’ 

Pag it) Pgg it) . 


(30) 


where analytic expressions of all the elements of the density matrix in Eq. (l30l > can be seen from Eqs. (23)-(32) of Ref. 172]. 

Here, we consider identical qubits. The dynamics involve collective coherent effects due to the multiqubit interaction, as 
well collective incoherent effects due to dissipative multiqubit interaction with the bath, and spontaneous emission. Analytic 
expressions of the corresponding QDs are very cumbersome, hence we resort to numerically studying the QDs for some param¬ 
eters. Values of different parameters are as follows wavevector and mean frequency ko = ojo = 1, spontaneous emission rate 
10 = 0.05, and p ■ fij = 0, where p is equal to the unit vector along the atomic transition dipole moment and f,y is the inter¬ 
atomic distance. Considering the initial state with p ee (0) = p gg (0) = p es (0) = p ea (0) = p eg (0) = p sg (0) = p ag (0) = 0, 
and p ss (0) = p aa (0) = p sa (0) = 0.5, the W, P, and Q functions are calculated. 

The variation of the different QDs is depicted in Figs. 0and0 Figs. 0a and 0 b show the negative values of W function and 
P function for some and for all times, respectively. The Q function exhibits a decaying pattern. These features are reinforced in 
the last plot of the figure, where all the QDs are plotted together. In Fig. 0 various QDs are plotted with respect to the inter-qubit 
distance. In the collective regime, r\ 2 -C 1, the QDs exhibit an oscillatory behavior, in consonance with the general behavior in 
this regime |t|]. Also, for the chosen parameters, the P function is always negative, while the W function is negative for t = 1, 
but becomes positive for a longer time t = 5, due to the dissipative influence of the bath. 


2. Squeezed thermal bath 

Fet us consider the evolution of same initial state, as in the presence of the vacuum bath at time t = 0, but now evolving 
under the influence of a squeezed thermal bath with finite T and r. Similar to the case of the vacuum bath, for certain values 
of the parameters in the density matrix, we can calculate different QDs. The elements of the density matrix in the presence of a 
squeezed thermal bath can be obtained as in Eqs. (33)-(40) in Ref. 117211 . For simplicity, we take here the squeezing angle $ = 0, 
and all other parameters are same as in the case of the vacuum bath, i.e., ko = cog = 1, Ty = 0.05, and p, ■ fij = 0. From 
the density matrix of the evolved state, the various QDs can be obtained, which once more due to their cumbersome nature are 
studied numerically. The behavior of the different QDs is shown, for different parameters, in Figs. 0and0 

In Figs. 0 three dimensional plots of W and P QDs are shown with respect to the azimuthal angles. The P function exhibits 
negative values for all values of the parameters chosen, while the W function does so for a restricted set of values. All these 
reiterate the quantumness of the state studied. From Fig. 0a-b, the effect of finite bath squeezing r and T on the evolution of 
the QDs can be seen. In particular, with an increase in T, the QDs, both P and IV', which were earlier exhibiting negative values 
start becoming positive, a clear indicator of a quantum to classical transition. Fig. 0c and d, showing the behavior of the QDs 
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Figure 4: (Color online) The variation of the W, P, and Q functions with time is shown in (a)-(c) for the two-qubit state, in the presence of 
vacuum bath, with 9i = 62 = </>i = j, <jf >2 = f with the inter-qubit spacing r 12 = 0.05 (smooth blue line), and ri 2 = 2.0 (red dashed 

line). In (d) all the QDs, varying with time, are plotted together with inter-qubit spacing ri 2 = 0.05. Here the P function is seen to be negative 
for all times shown, while the W function is negative only for sometime, while the Q function is always positive. 
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(b) 

QDs 



hi 




Figure 5: (Color online) (a) and (b) depict the W and P functions for the two qubit state, interacting with a vacuum bath, as a function of the 
inter-qubit spacing at t = 1 (smooth blue line), and t = 5 (red dashed line). In (c) and (d), W (smooth blue line), P (red dashed line) and Q 
(magenta dot-dashed line) QDs are plotted together, depicting their variation with inter-qubit spacing at time t = 1 and t = 5, respectively. 
For all the plots 61 = |, 62 = f, 4>i — f, and 4 >2 = f ■ 
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Figure 6 : (Color online) W and P functions were able to detect the nonclassicality in the presence of squeezed thermal bath at time t = 2.0 
with 9i = 62 = at temperature T = 1.0, squeezing parameter r = 0.5, and interatomic distance kr 12 = 0.05; all the remaining 

parameters are as mentioned above. The Q function failed to detect nonclassicality. The P function shows nonclassicality for all values of <j> 1 
and <j >2 as it is negative for all values of the chosen parameters. 


with respect to time and inter-qubit separations, are also along the expected lines. The nature of the QDs for small interqubit 
spacing, as seen here, is consistent with that of the same state subjected to dissipative interaction with a vacuum bath. Due to 
increase in temperature the oscillations observed in the QDs for small interqubit spacing, for the case of interaction with the 
vacuum bath in Fig. |5Jl, decreases in the present scenario of a squeezed thermal bath interaction and depicted in Fig. [7Jl. Further 
increase in temperature flattens the peak, observed for small inter-qubit spacing, in Fig. [7Jl. 


3. EPR singlet state in an amplitude damping (AD) channel 


Now, we take an initially entangled two qubit, EPR singlet, state |73i l. The evolution of this state is studied assuming inde¬ 
pendent action of an amplitude damping (AD) channel on each qubit. Such a scenario could be envisaged in a quantum memory 
net with the qubits being its remote components, subject locally to the AD noise (74]]. Using Kraus operators of an AD channel 


E 0 


V 1 - A(i) O' 

0 1 ’ 


Ei 


0 O' 

Vm oj ’ 


(31) 


where A (t) = 1 — e _7ot , where 70 is the spontaneous emission rate, and assuming that the two qubits, of the singlet, are 
independent and do not have any interaction, the Kraus operators for the action of AD channels, one on each spin, can be 
modeled as 


Ki = E 0 (A)®E 0 (B), 

Ix 2 = E 0 (A) ® Ei (B ), 

Ix 3 = Ei (A) <g> E 0 (B ), 

K 4 = Ei (A) ® Ei (B), 

where A and B stand for the first and second qubits (spins) comprising the singlet, respectively. From the form of Kraus 
operators (13T1 > and assuming A /( = A/> = A, we have 

' 1 - A 0 0 O' 

_ 0 vT - A 0 0 

Vl 0 0 sfpE A 0 ’ 

0 0 0 1 


0 0 0 O' 

sj\(\ - A) 0 0 0 

0 0 0 0 ’ 

0 0 Vx 0 
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Figure 7: (Color online) Different QDs, evolving under the influence of a squeezed thermal bath, are plotted with 0i = 62 = g, (f>i = 

f, <^2 = j, where smooth (blue), dashed (red) and dot-dashed (magenta) lines correspond to the W, P and Q functions, respectively, (a) and 
(b) are for temperature variation at time t = 2.0 in two cases, (a) r = 0.5, kr 12 = 0.08 (collective model), and (b) r = —0.5, kr\2 = 1.5 
(localized model). In (c) all the QDs are plotted with respect to time with T = 1.0, r = 0.1, kr\2 = 0.05. And (d) shows the variation of all 
the QDs withri 2 for T = 1.0, t = 2.0, r = 0.5. 


K-i 


0 0 0 O' 

0 0 0 0 

^(1 - A) 0 0 0’ 

0 v'AOO 


' 0000 ' 
0 0 0 0 
0 0 0 0 
.A 0 0 0. 


(32) 


The density matrix of the singlet state at time t, under the action of the above channel is 

pit) = 

1=1 


where p( 0 ) = \(f>) ((f >|, and 
density matrix is 


1 


■ 5 ) — | — 5 , 5 )) , is the initial state at time t = 0. Hence, at time t the evolved 


0 0 0 0 
0 |(1-A) -1(1 —A) 0 
0 -1(1-A) 1(1 -A) 0 
0 0 0 A 


(33) 


On the evolved state, represented by the above density matrix, we may now apply the prescription for obtaining the QDs to 
yield compact analytical expressions of the various QDs. Specifically, the W function is obtained as 

W (0i, 0i, 62 , (fe) = [A {1 + 3 cos 9\ cos 62 — \/3 (cos 9i + cos # 2 )} 

+ (1 — A) {1 — 3 cos 0i cos 02 — 3 sin0i sin0 2 cos (</> 1 — </> 2 )}], 


( 34 ) 
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Figure 8 : (Color online) The variation of all the QDs with time t for EPR singlet state in the presence of AD channel with 70 = 0.1 and 
8 \ = j, 82 = ■§, <t>i = 7,02 = Smooth (blue), dashed (red) and dot-dashed (magenta) lines are for W, P and Q functions, respectively. 

while the P function is 

P (0i, </>i, 0 2 , 02 ) = Yg^p- [A {1 + 9 cos 81 cos 02 + 3 (cos 6 \ + cos0 2 )} 

+ (1 — A) {1 — 9cos$i cos0 2 — 9sin0i sin0 2 cos (0i — 0 2 )}], ' 

and the Q function is 

<9(01,01,02)02) = [A{1 + COS01 COS02 + (COS01 + COS0 2 )} 

+ (1 — A) {1 — cos0! cos0 2 — sin0! sin0 2 cos (0! — 0 2 )}]. ' 

The QDs, reported here for an EPR pair (singlet state) evolving under AD channel exactly match with the corresponding noiseless 
results 11441 14911 . by setting A = 0 in the above expressions. The variation of the different QDs with time is shown in Fig. [ 8 ] 
The P and W functions are found to show negative values for a long time, indicative of the perfect initial entanglement in the 
system, and finally become positive due to exposure to noise. 

As all the QDs, in this case, are symmetric pairwise in (0i 7T 0 2 ) and (0 1 -f-f 0 2 ), hence either or both of these exchanges 
would leave the expressions unchanged. For 0\ = —0 2 = and 0i = 0 2 = 0, we observe the classically perfect anti¬ 
correlation of spins. We can also observe that for certain angles, viz. 0i = 0 2 = and Qi - 0 2 = where n is an odd 
integer, all the QDs become equal to , a result which remains unaffected by the presence or absence of noise. Hence, for 
these settings, the evolution of the QDs becomes noise independent. 


C. Three qubit QDs evolution in an AD channel 

Three qubit entangled states can be classified into two classes (GHZ and W classes) of quantum states, such that a state of W 
(GHZ) class cannot be transformed to a state of GHZ (W) class by using FOCC (local operation and classical communication). 
Here, we study both GHZ [75]] and W [76] classes of states. To simulate the effect of noise, we consider the scenario wherein 
the first qubit is affected by the AD channel. An arbitrary effect of noise on each subsystem could be thought of as more natural. 
The assumption of only one subsystem affected by amplitude damping noise is consistent with the effect of noise considered 
in, for example, various cryptrographic protocols ( fvh and references therein), where it is commonly assumed that the qubits 
which travel through the channel are affected by noise while the channel noise does not affect the qubits to be teleported or not 
travelling through it. 


1. GHZ state in an amplitude damping (AD) channel 

The GHZ (Greenberger-Horne-Zeilinger) state is a three qubit quantum state \GHZ) = (|000) + |111)). The first 

qubit of the state is acted upon by an amplitude damping (AD) channel while the remaining two qubits remain unaffected. The 
Kraus operators for AD channel for a single qubit state are as in Eq. (im Here, assuming that the three qubits are independent 
of each other and do not have any interactions, the Kraus operators for the action of AD channel only on the first qubit, can be 
modeled as 


Ki = E 0 (A) (g> I (B) (g> I (G), 
K 2 = Ei (A) ® I (B) <g) I (C), 


( 37 ) 
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where Eq and E\ are as in Eq. (OH and I is a 2 x 2 identity matrix. Also, A, B and C stand for the first, second and third qubits 
of the GHZ state, respectively. Thus, the density matrix of the GHZ state at time t in the amplitude damping channel is 


P (t) = hKi{t)p{Q)K\{t), 
2=1 


(38) 


where p (0) = \GHZ) ( GHZ \, is the initial state at time t = 0. Thus, at time t the density matrix for the GHZ state evolving in 
the presence of the AD channel is 


P(t) = 


(1-A) 

0 

0 

0 

0 

0 

0 


0 0 0 0 0 0 ,/(! - A) 


0 0 0 0 0 0 
0 0 0 0 0 0 
0 0 0 0 0 0 
0 0 0 A 0 0 
0 0 0 0 0 0 
0 0 0 0 0 0 


. sj(l - A) 0 0 0 0 0 0 

Analytical expressions can be obtained for the different QDs for the time evolved GHZ state described by Eq. d39l ). 
Specifically, we obtain the W function as 

W ( 01 , 01 , 02 , 02, $ 3 , 03) = 0 ^ 3 - [l — -\/3AcOS0i + 3 COS 02 COS $3 + 3 (1 — A) COS 01 (cos 02 + COS 03 ) 


0 'j -L y U / \ vv V_Z U l J y \ Kj v/WO 1/ ^ Vv V/ U 1/ ^ | KJ y JL / \) v/C/U 1/ ^ 1/ ^ | 

3x/3 j A cos 0i cos 02 cos 0 3 — ^/(l — A) sin 0i sin 0 2 sin 0 3 cos (0 1 


■ 0 2 


while the P function is 


P (01, 01,02, 02 , 03 , 03) = [1 + 3ACOS01 + 9 COS 02 COS 03 + 9 (1 - A) COS01 (cos 02 + COS0 3 ) 

■ 02 + 


647T 

+ 27 


| A cos 0i cos 02 cos 0 3 — i/(l — A) sin 0i sin 0 2 sin 0 3 cos (0i 


and the Q function is 


Q (01,01,02, <(>2,03, 0 3 ) = 647^ I 1 + A COS 01 + COS 02 COS 0 3 + (1 - A) COS 01 (COS 0 2 + COS0 3 ) 


+ A cos 0i cos 02 cos 0 3 — -\/(l — A) sin 0i sin 02 sin 0 3 cos (0i + 02 + 0 3 ) 


(39) 


(40) 


(41) 


(42) 


The variation of the different QDs, as in Eqs. (|401>- (|4H . with time is depicted, for a particular choice of the parameters, in Fig. 
[9}t. The W and P functions exhibit negative values (nonclassical character) for the times shown, which could be attributed to 
the initial entanglement in the state. Further, for the sake of generality, we depict in Fig. |9j> the scenario wherein all three 
qubits are affected by the generalized amplitude damping noise lf58|] , which can be obtained from Eq. (l24l > by setting the bath 
squeezing parameters to zero. All three qubits are subjected to different temperatures corresponding to independent environment 
for each qubit, and simulates the scenario where each qubit travels through an independent channel. It can be observed that the 
nonclassicality indicated by the negativity of W and P functions at t = 0 decays more rapidly when the last two qubits are 
subjected to finite temperature noises. 


2. W state in an amplitude damping (AD) channel 


For our second example of QDs of three qubit states, we take up the W state | W) — (1001) + |010) + 1100)). As 

before, we consider the evolution where only the first qubit of the state is acted upon by an AD channel. The Kraus operators, 
describing the evolution are given by Eq. (1371) . Here, as in the case of GHZ or EPR states, assuming that the three qubits are 
independent of each other and do not have any interactions, the density matrix of the evolved state is, as in the last case, given 
by Eq. ( |38] >. where p (0) = | W) {W \, is the initial state at time t = 0. Hence, at time t, the W state evolves, in the presence of 
the AD channel, to 


P(t) 


0 0 


0 0 0 0 0 0 


(1-A) 

(1-A) 

0 


(1-A) 0 a/(1 - A) 0 0 0 
(1-A) 0 - A) 0 0 0 

0 0 0 0 0 0 


o a/0~a) vQ^a) 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 


1 

0 

0 

0 


0 0 0 
A A 0 
A A 0 
0 0 0 


1 _ 

3 


(43) 
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Figure 9: (Color online) The variation of all QDs with time for the GHZ state when acted upon by (a) an AD noise on the first qubit and (b) 
a generalized amplitude damping noise on each qubit, with 70 = 0.1 and 9\ = ^,82 = §, 83 = j, <fii = §2 = f 1 03 = f ■ In (b) 
different QDs are shown for w = 1.0 with the first, second and third qubits subjected to generalized amplitude damping noise at T = 0, 1 and 
2, respectively. In both the plots smooth (blue), dashed (red) and dot-dashed (magenta) lines are for the W, P and Q functions, respectively. 


In this case, again, making use of Eq. (I43l) . analytical forms of the different QDs can be obtained as follow; 

W (01,01, 02, 02,03, 0 * 3 ) = 1 - (cos di COS d 2 + COS 0 2 cos 0 3 + COS 9i COS 0 3 ) + ^ (cos 9i + 

~ “ " - 2 (l +y/3cos0i) sin 02 sin 0 3 cos (02 — 0 3 ) 

3 82 ) sin 8 1 sin 0 3 cos (0i — 0 3 ) 


3 -S/ 3 COS 01 cos 6*2 cos 0 3 z (i - 1 - vocostfij sin 02 sine's cos {(P 2 — 
2a /(1 — A) { (l + y/3c O s0 2 ) sin sin0 3 cos (0i — 0 3 ) 

(l + -\/3 cos 0 3 ) sin 0 i sin 0 2 cos (0i — 02)} 

4\/3Acos0i { — | + cos 02 cos0 3 — sin0 2 sin0 3 cos (02 — 03)}] , 


and 


P ( 01 , 01 ,02 , 02, 03, 03) = 0^3 [1 — 3 (cos 01 COS 0 2 + COS 0 2 COS 0 3 + COS 0i COS0 3 ) — (COS 0i + COS 0 2 
+ 27 cos 0i cos 02 cos 0 3 + 6 (1 — 3cos0i) sin0 2 sin0 3 cos (02 — 0 3 ) 

+ 6y/(1 — A) {(1 — 3 cos02) sin 0i sin0 3 cos (0i — 03) 

+ (1 — 3cos 0 3 ) sin 0 i sin 0 2 cos (0i — 02 )} 

+ 4A cos0i {1 — 9 cos02 cos 0 3 + 9 sin0 2 sin0 3 cos (02 — 03)}], 


- cos 0 3 ) 


Q ($ 1 ) 01, 82 , 02 


, 03 , 03) = 19 } 7r 3 [3 ~ (COS 0i COS 02 + COS 02 COS 0 3 + COS 01 COS 0 3 ) — (COS 01 + COS 02 

Q 2 

+ 3 cos 0i cos 0 2 cos 0 3 + 4 sin -f sin 0 2 sin 0 3 cos (02 — 03) 

+ 4yj (1 — A) < sin0i sin sin0 3 cos (0i — 0 3 ) + sin0i sin0 2 sin cos (0; 
+ 4A cos 0 i {1 — cos 02 cos 63 + sin 0 2 sin 0 3 cos (0 2 — 03 )}] ■ 


cos 83 ) 

- h)} 


(44) 


(45) 


(46) 


The variation of the QDs with time is shown in Fig. ITOh for a particular choice of the parameters. Here, the P function exhibits 
negative values for the times shown, but in contrast to the GHZ case, the W function is found to be positive. Thus, the signature 
of quantumness (nonclassicality) of the state is identified by P function, but W function fails to detect the same. Similar to the 
QDs of GHZ state in the presence of generalized amplitude damping noise, a decrease in the nonclassicality of W state with 
time can be observed due to the effect of noise on the last two qubits at non-zero temperatures (cf. Fig. ITOb). 


D. N qubit Dicke Model 

We conclude our discussion of QDs for two-level systems (qubits) with the Dicke model. The interaction of a collection of 
identical two-level atoms with a single mode of quantized field is the Dicke model ll22U78ll and is the multi-atom generalization 
of the Jaynes-Cummings model. The Hamiltonian for the Dicke model is given as 

H = w(a t a + E°) + <7 (a£^ + a^Ei) . (47) 

Here, ui is the resonant atomic frequency, g is the coupling constant, a («') is the annihilation (creation) operator for the 

N N 

radiation field, and E[j_ = and E( = where N is the number of atoms considered. Also, the symbol c 

3 =1 i=i 
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Figure 10: (Color online) (a) Variation of QDs with time for the W state in the presence of an AD channel, acting only upon the first qubit, 
with 70 = 0.1 and 6\ = 62 = 63 = (j>\ = j, 4 > 2 = f, <t>3 = • (b) QDs for W state, when all the qubits are subjected to 

generalized amplitude damping noise, with uj = 1.0 and T = 0, 1 and 2 for the first, second and third qubits, respectively. The remaining 
values of parameters in (b) are as in (a). Smooth (blue), dashed (red) and dot-dashed (magenta) lines correspond to W function, P function 
and Q function, respectively. 


QDs 



Figure 11: (Color online) Variation of all QDs with time for the IV = 4 atoms Dicke model with 6 = and cj> = Smooth (blue), large 
dashed (red), small dashed (magenta) and dot-dashed (cyan) lines are for W function, P function, Q function, and F function, respectively. 
The four atom case is equivalent to the spin-2 case and consequently different values of the W and F functions can be observed. 


in the superscript stands for collective. A generalization of this model, dynamics of a collection of atoms interacting with a 
squeezed radiation field, was made in |79]. The Dicke model has analytical solutions in two regimes: weak and strong field 
regimes, which corresponds to average photon number being much smaller or greater than the number of atoms, respectively. 

We discuss, here, only the second case in the presence of strong initial field l80ll . i.e., ii> N, where n is the average number 
of photons in the initial coherent field. In the dissipative case, modeling the microwave region of zero temperature dissipative 
cavity quantum electrodynamics, restricting ourselves to the condition ^ <C \pfi, the master equation in Fock basis is 

Pnm 21(7 Pnm Pnm N £) T ^ (2 ClpO(L* dp p(l^ , (48) 

where un = n + \ and 7 is the cavity decay constant. This master equation takes a simple form in the dressed atomic 
basis, in which the atomic operator \ (£(j_ + £i) is diagonal 

££ \q) = K \q), 

where X q = q -j ■ an d 9 = 0,1,... ,N. These dressed atomic states can be expressed in terms of the bare atomic basis as 
l<?> = J2Cqk |fc), where C qk = (k |g) = i q ~ k d* k (-§) , and d qk (6) are Wigner d functions. The bare atomic basis (Dicke 

Q 

states) is defined as 

E'l*> = (*-y) l fc >- 0 <k<N, 


where k is the number of excited atoms. 
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Now, we consider the initial atoms-field density matrix 

p{ 0 ) = |atin) (atin| ® |a) (a|, 

where |atin) is the initial atomic state, taken here to be ground state, i.e., none of the atoms are excited and a) denotes the initial 
strong coherent state of the field. The atomic density matrix, in the bare atomic basis, after tracing out the radiation field evolves 
as 


Pki (t) 


^jOik+n (t) a l-\- n 

n 


N 


(t) E C kq C lp C 0 q C* 0 p G%+ k ’ n + l (t), 

q,p —0 


(49) 


where 


G qp m ( t) = exp [-2 igt (A qy /n^ - A py /m^) - O qp (t)], 

and O qp (t) = n |(1 - e -7 *) - (l - e -7 ' 4 ) j, V = 7 + and “» {t) = e ~ ra(4)/2 , with n ( t) = ne _7t . 

The evolution of all the QDs can be obtained for the Dicke model in the dissipative case using Eq. ( l49l >. Here, we consider 
the N = 4 atoms Dicke model in a strong initial field with average photon number 30. Further, we take A = 10” 2 , and coupling 
constant as 0.1. To achieve consistence with the notations used in this article, we have used j = -j, and m, m' = —j,... , j. 
The multipole operators required to calculate the QDs are equivalent to spin-2 multipole operators and can be seen in Appendix. 
In Fig. QT| variation of the various QDs, as they evolve with time, is depicted. The P function exhibits negative values at some 
early times, indicative of the quintessence of quantumness in the system, but eventually becomes positive due to dissipative 
effects. Different values of the W and F QDs are observed at some time intervals, consistent with the observation that the four 
atom case is equivalent to the spin-2 case. For this specific choice of parameters and restricting ourselves to the computational 
accuracy of the numerical method adopted here, we did not observe any signature of nonclassicality via W function. However, 
W function can witness the nonclassical characteristics present in the Dicke model for other values of 0 and </>. This point will 
be clearly illustrated in Fig. [13]d, where we show positive values of nonclassical volume (which is obtained by integrating the 
modulus of W function over all possible values of 0 and 0 ) for the Dicke model. These positive values of nonclassical volume 
imply negative values of W function for some values of 9 and 0 (cf. (Eg. [54b. 


V. QDS FOR A SPIN-1 STATE 

i from spin-A to spin-1 states. For a : 

lx) = °+ 1+1) + °o |0> + a_ |—1), 


Now, we extend the discussion of spin QDs from spin-1 to spin-1 states. For a spin-1 pure state 


T °° - 73 
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. All other multipole operators can be obtained from these operators. The analytical expressions of the 


different QDs are obtained as 

W(9,(/)) = y+jt 4 — v/l0 + 3-\/l0 ^cos 2 9 + |ao | 2 — 3 |ao | 2 cos 2 9^j + 6v/2 

+ ( 6 ao sin 0 {a?j_ exp(— i(j>) (l + v/Scosf?) + a*_ exp(l — \/5cc 


|a +| 2 — |a_| 2 ) cost 


cos 9) } + 3VT0a+ail sin 2 9 exp( 2 i(f>) + c.c.)] , 

(50)' 

. cos 9 

0)} + 5a+al sin 2 0exp(2*^>) + c.c.)] , 

<: 


P{0,4>) = -1 + 5 (^cos 2 9 + |ao | 2 - 3 |a 0 | 2 cos 2 9^j - 2 (|a +| 2 - |a_| 2 ^) 

+ (v/2ao sin 0 {a+ exp(—*</>) (1 — 5 cos 0) + a* exp(i^) (1 + 5 cos t 

Q{9,cj)) = 1 + ^cos 2 0 + |ao|“ — 3 |ao | 2 cos 2 0^ — 2 ^|a +| 2 — |a_| 2 ^ cos0 

+ (v/2ao sin0 {a^ exp(— i(j>) (1 — cos0) + a* exp(i<^) (1 + cos0)} + 5a + al sin 2 0exp(2i<^) + c.c.)] , 


(51) 


( 52 ) 
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QDs QDs 




Figure 12: (Color online) Variation of all the QDs is shown for the spin-1 state with a+ = ao = a_ = ^=. (a) shows W (smooth blue line), 
P (small dashed red line), Q (large dashed green line), and F (dot-dashed magenta line) functions with <j> = Similarly, (b) shows all the 
QDs for 9 = ^. Both (a) and (b) clearly bring out the point that W and F functions are not equal, in general, in contrast to the spin-4 cases. 


and 


F(M) = 
+ 


gf^r 1 + 5 ^cos 2 9 + |ao| 2 — 3 |ao| 2 cos 2 9^j + Ay/2 ^|a+| 2 — |a_| 2 ^ cos# 

(ao sin# {a?J_ exp {—i<f>) (4 + 5-\/2cos 9) + a* exp(4 — 5-\/2cos 9 )} + 5a + al sin 2 #exp(2i</>) + c.c.)] . 


(53) 

In Fig. [T21 we illustrate the behavior of the above QDs for a_|_ = ao = a_ = -4=. Here, we do not consider the effect of noise on 
the evolution of the QDs, a topic to which we will return back to, in the future. The purpose here, besides studying the quantum¬ 
ness in the system via the QDs, is to emphasize the nonequivalence, for the spin-1 case, of the W and F functions, in contrast 
to the spin-4 case. Fig. [I2]depicts the behavior of the QDs with respect to 9 and <j>, both of which show a symmetric behavior 
about the central point on the ordinate. The P, W and F functions exhibit negative values, indicative of the quantumness in the 
system, with the P function being the most sensitive indicator, as expected. Also, the F and W functions are clearly distinct. 


VI. NONCLASSICAL VOLUME 

Till now, we have studied nonclassicality using negative values of the W or P function. Negative values of the QDs only 
provide a signature of nonclassicality, but they do not provide a quantitative measure of nonclassicality. There do exist some 
quantitative measures of nonclassicality, see for example, lf68tl for a review. One such measure is nonclassical volume introduced 
in [69]. In this approach, the doubled volume of the integrated negative part of the W function of a given quantum state is used 
as a quantitative measure of the quantumness [69|]. Using our knowledge of the W functions for various systems, studied here, 
the nonclassical volume 6 , which is defined as 


S = J \W (9, (j>)\sm9d9dcp — 1, (54) 

can be computed. It can be easily observed that a nonzero value of 5 would imply the existence of nonclassicality, but this 
measure is not useful in measuring inherent nonclassicality in all quantum states. This is so because, the W function is only a 
witness of nonclassicality (it does not provide a necessary condition). However, this measure of nonclassicality has been used in 
a number of optical systems, see for example, l(8lll82ll and references therein. 

Here, we will illustrate the time evolution of 5 for some of the spin-qubit systems, studied above. Specifically, Fig. Q2]a 
and b show the variation of S for two spin- 4 systems, initially in an atomic coherent state under the influence of QND and 
SGAD channels, respectively. Fig. [I3]c shows two spin-4 states in a two-qubit vacuum bath. The dashed and dot-dashed lines 
in Fig. [13] a and b, i.e., for atomic coherent state in QND with finite temperature and in SGAD with finite temperature and 
squeezing, exhibit the exponential reduction of nonclassical volume with time implying a quick transition from nonclassical to 
classical states, whereas the smooth lines in a and b (when temperature and squeezing parameters are taken to be zero) show 
that after an initial reduction, the nonclassical volume stabilizes over a reasonably large duration. Thus, nonclassicality does not 
get completely destroyed with time. A similar nature of time evolution of 6 is also observed for the dashed line in Fig. [T3]c (a 
two qubit state in a vacuum bath with relatively large inter-qubit spacing), whereas an oscillatory nature is observed for small 
inter-qubit spacing, depicted here by a smooth line. It should be noted that the smooth blue line in Fig. Qj] b corresponds to 
the nonclassical volume for an atomic coherent state dissipatively interacting with a vacuum bath, i.e., at zero temperature and 
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5 6 



Figure 13: (Color online) The plots (a)-(d) depict the variation of the nonclassical volume in the presence of different noises, (a) The 
variation of the nonclassical volume with time is shown for the single spin-^ atomic coherent state in the presence of the QND noise with 
70 = 0.1, r = 0, a = 0, ui c = 100, u> = 1.0, and a = j, (3 = 7 , in the units of h = Icb = 1, where the smooth (blue), dashed (red) and 
dot-dashed (magenta) lines correspond to different temperatures T = 0, 1 and 2, respectively, (b) The variation of the nonclassical volume 
with time is shown for a single spin-i atomic coherent state in a SGAD channel, where the smooth (blue) line corresponds to the variation 
in nonclassical volume in vacuum bath, i.e., at T = 0 and squeezing parameters r = £ = 0 (amplitude damping channel); dashed (red) line 
corresponds to the variation in a channel with zero squeezing at T = 3, i.e., generalized amplitude damping channel, and dot-dashed (magenta) 
line corresponds to the variation with squeezing r — 1, squeezing angle £ = 0, and T = 3. In all these cases a and f3 have the same values 
as in (a), (c) The behavior of nonclassical volume with time is depicted in a vacuum bath for the state discussed in Section ITV B II with the 
inter-qubit spacing r 12 = 0.05 (smooth blue line), and 772 = 2.0 (red dashed line), (d) The temporal behavior of nonclassical volume of the 
four atom Dicke model interacting with a strong input coherent field with average photon number 30, in the dissipative regime with ^ = 10~ 2 , 
and coupling constant as 0 . 1 . 


squeezing, while the W function of the atomic coherent state, illustrated in Fig. [2] is for non-zero temperature and squeezing. At 
zero temperature, the nonclassicality present in the system is expected to survive for a relatively longer period of time. Further, 
the nonclassical volume is the overall contribution in nonclassicality from all values of 6 and <f>. It is possible that the W function 
shown for a particular value of 9 and </>, in the previous sections, may not exibit nonclassical behaviour at time t whereas the 
other possible values provide a finite contribution to nonclassical volume, resulting in a nonvanishing 5. Interestingly, in Fig. [13] 
d for the nonclassical volume of the Dicke model, we find that the amount of nonclassicality oscillates with time. This seems to 
arise from the weakness of the nonclassicality measure used here. 


VII. CONCLUSIONS 

The nonclassical nature of all the systems studied here, of relevance to the fields of quantum optics and information, is illus¬ 
trated via their quasiprobability distributions as a function of the time of evolution as well as various state or bath parameters. We 
also provide a quantitative idea of the amount of nonclassicality observed in some of the systems studied using a measure which 
essentially makes use of the W function. These issues assume significance in questions related to quantum state engineering, 
where the central point is to have a clear understanding of coherences in the quantum mechanical system being used. Thus, it 
is essential to have an understanding over quantum to classical transitions, under ambient conditions. This is made possible by 
the present work, where a comprehensive analysis of QDs for spin-qubit systems is made under general open system effects, 
including both pure dephasing as well as dissipation, making it relevant from the perspective of experimental implementation. 
Along with the well known W, P and Q quasiprobability distributions, we also discuss the so called F function and specify 
its relation to the W function. We expect this work to have an impact on issues related to state reconstruction, in the presence 
of decoherence and dissipation. These quasiprobability distributions also play an important role in fundamental issues such as 
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complementarity between number and phase distributions as well as for phase dispersion in atomic systems. It is interesting to 
note that in Issll . a connection was established between negative values of a particular quasiprobability and potential for quantum 
speed-up. The present study could be of use to probe this connection deeper. 
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Appendix: Multipole operators for Dicke model 


Here, we collect the multipole operators used for the computation of QDs for the four qubit Dicke model in Sec. (IV.D). 
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